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Bats have been implicated as important reservoir hosts of alpha- and betacoronaviruses. In this 
study, diverse coronaviruses (CoVs) were detected in 50 of 951 (positive rate 5.3%) intestinal 
specimens of eight bat species collected in four provinces and the Tibet Autonomous Region of 
China by pan-coronavirus RT-PCR screening. Based on 400-nt RNA-dependent RNA polymerase 
(RdRP) sequence analysis, eight belonged to genus Alphacoronavirus and 42 to Betacoronavirus . 
Among the 50 positive specimens, thirteen gave rise to CoV full-length RdRP gene amplification for 
further sequence comparison, of which three divergent sequences (two from a unreported 
province) were subjected to full genome sequencing. Two complete genomes of betacoronaviruses 
(JTMC15 and JPDB144) and one nearly-complete genome of alphacoronavirus (JTAC2) were 
sequenced and their genomic organization predicted. The present study has identified additional 
numbers of genetically diverse bat-borne coronaviruses with a wide distribution in China. Two new 
species of bat CoV, identified through sequence comparison and phylogenetic analysis, are proposed. 
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INTRODUCTION 

Coronaviruses (CoVs, family Coronaviridae , subfamily 
Coronavirinae) are important human and animal patho¬ 
gens which, according to the latest release of Virus Tax¬ 
onomy by the International Committee on Taxonomy of 
Viruses (ICTV, http://www.ictvonline.org/virusTax- 
onomy.asp?msl_id=26), currently comprise four distinct 
genera: Alphacoronavirus ( aCoV ), Betacoronavirus 
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( jSCoV ), Gammacoronavirus ( yCoV) and Deltacronavir- 
us (SCoV ). This large group of viruses has a wide spec¬ 
trum of hosts, including humans, rodents, carnivores, 
chiropters and avians, and cause respiratory, enteric, hep¬ 
atic and neurological diseases (Lai et al, 2007). They in¬ 
clude even public threats such as the severe acute respir¬ 
atory syndrome (SARS) and the current Middle East res¬ 
piratory syndrome (MERS) (Moratelli et al, 2015). Bats 
are host animals of diverse aCoVs and pCoVs that may 
serve as the ancestral origins of mammalian Co Vs (Fal¬ 
con et al, 2011; Woo et al, 2012). In last decade, increas¬ 
ing numbers of bat Co Vs with wide molecular diversit¬ 
ies have been reported worldwide, particularly in China 
(Li et al, 2005; Tang et al, 2006; Woo et al, 2007; Chu et 
al, 2008; Yuan et al, 2010; He et al, 2014), some of 
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which likely have the potential ability to cause human 
diseases (Ge et al, 2013; He et al, 2014; Menachery et al, 
2015). These findings indicate that further diverse CoVs 
circulate in bat populations. China has a nationwide dis¬ 
tribution of about 120 bat species, with many roosting re¬ 
gions remaining uninvestigated for harbored mammal 
viruses. Here, we report a continuing investigation on 
bat-borne Co Vs in some unexplored regions in China, 
the results of which have revealed more novel Co Vs that 
circulate and evolve in bat populations with great mo¬ 
lecular diversity and wide geographic distribution. 

MATERIALS AND METHODS 

Bat collection and species conformation 

A total of 951 bats covering 5 families and 21 species 
were captured between 2005 and 2013, in Jilin, Liaoning, 
Yunnan, Guangdong provinces and the Tibet Autonom¬ 
ous Region, China. Bat species were morphologically 
identified by a trained field biologist and further con¬ 
firmed by PCR of their mitochondrial cytochrome b gene 
sequence (Wang et al, 2003). Respiratory and intestinal 
tissue specimens were collected separately from each bat 
and stored at -80 °C immediately until further pro¬ 
cessing. 

RNA extraction and detection by RT-PCR 

Viral RNA of each specimen was extracted by using the 
RNeasy Mini Kit (Qiagen, Hilden, Germany), and was 
immediately reverse-transcribed with the Superscript III 
Kit (Invitrogen, San Diego, CA) using random primers. 
Pan-CoV nested PCR primers were used to amplify a 
440-nt sequence in the RNA-dependent RNA poly¬ 
merase (RdRP) gene by our published methods (He et al. 
2014) (see Supplementary Table SI for primer informa¬ 
tion). Expected PCR amplicons were directly sequenced 
by the Sanger method in an ABI 3730 sequencer (Co¬ 
mate Bio, Changchun, China). 

RdRP gene amplification and whole genome 
sequencing 

The complete RdRP genes of positive samples were 
amplified using LA Taq (TaKaRa, Dalian, China). 
Primers were designed based on RdRP gene sequences of 
representative aCoV and fiCoV strains available in Gen- 
Bank. Reactions were carried out with a touch-down 
PCR program: 94 °C for 3 min, then 10 temperature 
decrement cycles (94 °C for 30 s, 58 °C minus 1 °C per 
cycle for 30 s, 72 °C for 2 min), followed by 35 normal 
cycles (94 °C for 30 s, 52 °C for 30 s, 72 °C for 3 min), 
and a last extension of 72 °C for 10 min. 

To obtain the full genomes of the interesting speci¬ 
mens, overlapping amplicons were obtained by the above 
PCR program following by assembly into contigs. In ad¬ 


dition, deep sequencing and genome walking were also 
undertaken to recover more genomic sequences. The 5' 
and 3' termini were sequenced using a 5' Full RACE Kit 
with TAP and a 3' Full RACE Core Set with PrimeScript 
RTase (TaKaRa, Dalian, China). Primer sequences for 
full-length genome amplification are shown in Supple¬ 
mentary Table SI. 

Genomic and phylogenetic analyses 

Genomic structures of the CoV complete sequences were 
predicted by the SeqBuilder program of the DNAStar 
software package and compared with other representat¬ 
ives from GenBank. Nonstructural proteins (nsps) in 
ORFla and ORFlb (replicase) of the Co Vs were pre¬ 
dicted using Z-Curve version 2.0, a CoV-specific gene¬ 
finding system (Gao et al, 2003). 

All 400-bp amplicons (the primer truncation of 440-nt 
sequences) were aligned with their closest phylogenetic 
neighbors in GenBank using Clustal W version 2.0. The 
phylogenetic tree was then constructed by the maximum 
likelihood method of MEGA 6.06 with 1,000 bootstrap 
replications. To better understand their evolutionary rela¬ 
tionships, the complete RdRP genes were further ampli¬ 
fied and used for the analysis. 

Nucleotide sequence accession numbers 

The partial RdRP sequences obtained from all positive 
samples and the complete genome or full length RdRP 
sequences of some specimens were submitted to the 
GenBank under accession numbers KU182954 to 
KU183005. 

RESULTS 

Detection of CoVs 

Of 951 bats tested 50 intestinal specimens (5.3%) were 
CoV positive, but surprisingly all respiratory specimens 
showed negative amplification. As shown in Table 1, 
among 181 bats from 6 species in 3 families in Guang¬ 
dong province, 16.2% (6/37) Rousettus leschenaultia and 
27.5% (14/51) Cynopterus sphinx were CoV positive. 
Among 599 bats from 17 species in 5 families in Yun¬ 
nan province, 14.0% (14/100) Rousettus leschenaulti , 
2.4% (1/41) Megaerops kusnotei , 9.0% (7/78) Rhinolo- 
phus sinicus and 5.3% (5/95) Myotis daubentonii were 
CoV positive. As the first study of this kind in the Tibet 
Autonomous Region, fifteen Hipposideros cineraceus 
and five Rhinolophus hipposideros collected in south 
Tibet were tested and only 6.7% (1/15) Hipposideros 
cineraceus showed positive amplification. In northeast 
China, 2 of 97 (2.1%) bats in Jilin province were posit¬ 
ive: one from Murina leucogaster and another from 
Rhinolophus ferrumequinum. In contrast, all 16 Rhinolo¬ 
phus ferrumequinum and 38 Myotis ricketti in Liaoning 
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Table 1. Bat sample collection and coronavirus detection. 


Bats 


Guangdong 
(Year 2005) 

Yunnan 

(Years 2012, 2013) 

Tibet 

(Year 2013) 

Liaoning 
(Year 2013) 

Jilin 

(Year 2013) 

Family 

Species 

Bat § 

CoV & 

Bat 

CoV 

Bat CoV 

Bat CoV 

Bat CoV 

Pteropodidae 

Rousettus 

leschenaulti 

6/37(16.2) 

(34 

14/100(14.0) 

(34 





Cynopterus 

sphinx 

14/51(27.5) 

(34 







Megaerops 

kusnotei 



1/41(2.4) 

(34 




Hipposideridae 

Hipposideros 

cineraceus 

0/9 




1/15(6.7) a 




Hipposideros 

pomona 



0/84 






Hipposideros 

larvatus 

0/68 


0/2 






Hipposideros 

armiger 

0/11 


0/18 






Aselliscus 

stoliczkanus 



0/33 





Rhinolophidae 

Rhinolophus 

ferrumequinum 



0/42 



0/16 

1/30(3.3) (32 


Rhinolophus 

sinicus 



7/78(9.0) 

a, (32 





Rhinolophus 

pusillus 

0/5 


0/6 






Rhinolophus 

affinis 



0/3 






Rhinolophus 

hipposideros 



0/37 


0/5 




Rhinolophus 

macrotis 








Vespertilionidae Myotis 

daubentonii 



5/95(5.3) 

(33 





Myotis laniger 



0/8 






Myotis 

chinensis 



0/3 






Myotis 

capaccinii 



0/40 






Myotis ricketti 






0/38 

0/27 


Miniopterus 

schreibersi 



0/8 






Murina 

leucogaster 







1/40(2.5) a 

Megadermatidae 

Megaderma 

lyra 



0/1 






p & 

Note: positive/total bats; numbers in brackets indicate the coronavirus positive percentage. CoV, a: aCoV; (3: unclassified 
(3CoV; (32: (3CoV lineage 2; (33: (3CoV lineage 3; (34: (3CoV lineage 4. 


province showed negative amplification. These results 
revealed a higher CoV incidence in three fruit bat spe¬ 
cies of the family Pteropodidae than in the four insecti¬ 
vorous bat families, indicating that fruit bats are more 
likely to harbor CoVs. 

Phylogenetic analysis 

To describe the genetic relationships among the 50 se¬ 
quences obtained in this study and previously known 


CoVs, 400-nt RdRP sequences were obtained from the 
primer truncation of 440-nt sequences and phylogenetic- 
ally analyzed. Results showed that 8 sequences grouped 
into 3 clusters within the genus aCoV (Figure 1A). YDB5C 
is the first reported bat-borne CoV (Hipposideros ciner- 
aceus ) in Tibet and clustered closely with MLHJC4, a 
CoV from Rhinolophus sinicus in Yunnan, both sharing 
94% nt identity with previously reported strain HKU2/ 
GD/430/2006 from Guangdong (Lau et al, 2007). JTAC2 
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identified in Murina leucogaster in Jilin province di¬ 
verged considerably from known Co Vs, showing the hig¬ 
hest nt identity of only 83% to bat-borne coronavirus 
Neixiang-14 and Neixiang-52 detected also in Murina 
leucogaster , and followed by 78% nt identity with some 
pandemic porcine epidemic diarrhea virus (PEDV) strains 
that have emerged recently in China, USA and Japan 
(Vlasova et al, 2014; Sun et al, 2015; Suzuki et al, 2015). 
Five other aCoVs (MLHJC1, MLHJC6, MLHJC8, ML- 
HJC22, MLHJC34) identified from Rhinolophus sinicus 
in Yunnan formed a new group with MLHJC8 being 
slightly more divergent, showing highest nt identities 
(75%-89%) with the previously reported BtCoV/860/2005 
(Tang et al, 2006). The remaining 42 bat CoV sequences 
were classified as fiCoV and fell into 5 clusters ( figure 
IB). Twenty identified in Guangdong fell into lineage 
P4, which showed the geographical relationship and was 
further divided into two distinct clusters, one with 6 se¬ 
quences sharing 99% highest nt identity with HKU9-10- 
1 (Lau et al, 2010), while another including 14 se¬ 
quences sharing the closest relationship with BtCoV/ 
BRT55630/H.lek/CK/Tha/05/2012 detected in Hippos- 
ideros lekaguli in Thailand (Wacharapluesadee et al, 
2015). The 21 pCoVs identified in Yunnan province ex¬ 
hibited considerable genetic diversity and were distrib¬ 
uted among lineages P2, P3 and p4. Fifteen fell into p4 
and further divided into 2 lineages, fourteen sequences 
showing closest relationship to previously reported Bt- 
CoV/BRT55629/H.lek/CK/Tha/05/2012 (Wacharapluesadee 
et al, 2015), while another (MF92C) grouped with the 
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Figure 1. The 400-nt RdRP gene fragment based 
phylogenetic analysis of 50 bat CoV sequences ob¬ 
tained in this study (8 sequences of aCoV in (A) and 42 
sequences of (3CoV in (B)) in comparison with other 
representative strains retrieved from GenBank. Fifty se¬ 
quences of this study are marked by triangles (13 se¬ 
quences with complete RdRP gene sequencing are 
marked by solid triangles). The scale bar indicates the 
estimated number of substitutions per 10 nucleotides. 
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above Guangdong sequences. Five sequences detected 
from Myotis daubentonii clustered within lineage P3, 
sharing > 91% nt identities with previously reported 
HKU4-4 from Tylonycteris pachypus (Woo et al, 2007). 
This group showed about 80% nt identity with MERS- 
CoVs recently identified in China (Lu et al, 2015) and 
Korea (Kim et al, 2015). The remaining Yunnan bat CoV 
sequence, MLHJC35, detected in Rhinolophus sinicus, 
and the only Jilin province sequence, JTMC15, identi¬ 
fied in Rhinolophus ferrumequinum were clustered into 
P2 and showed highest nt identities to SARS-related bat- 
borne CoVs (SARSr-BatCoVs). MLHJC35 was 97% 
identical with SARSr-BatCoV Cp/Yunnan2011 previ¬ 
ously isolated in Yunnan province (Yang et al, 2011), 
while JTMC15 shared 99% identity with SARSr-Bat- 
CoV Rfl found in Rhinolophus ferrumequinum, Hubei 
province (Li et al, 2005). 

To obtain more precise analysis, representative speci¬ 
mens of the 8 phylogenetic clusters were subjected to 
full RdRP gene amplification. Complete RdRP se¬ 
quences were obtained with 13 specimens belonging to 6 
clusters comprised of 3 aCoVs and 10 pCoVs. Phylogen¬ 
etic analysis based on the full RdRP gene sequences was 
highly consistent with figure 1 (Phylogenetic tree of the 
full RdRP gene sequences is not shown). 

Full genomic sequences characterization 

Full genomic sequencing was successful in 2 of the 
above 13 specimens: JTMC15 from Rhinolophus fer¬ 
rumequinum, Jilin province, and JPDB144 from My otis 
daubentonii, Yunnan province, with a nearly complete 
genome sequence obtained of JTAC2 from Murina leu- 
cogas ter. The full genomes of JTMC15 and JPDB144 
(including complete terminal sequences of 5' end and 3'- 
poly A) and near-complete genome of JTAC2 were 
28,761 nt, 30,321 nt and 25,719 nt in size respectively, 
with G+C contents of 38.1%, 41.0% and 43.4%. It is 
proved that two proteinases, papain-like proteinase (PL Pl °) 
encoded by nsp3 gene and main proteinase (M ) en¬ 
coded by nsp5 gene in ORFla of Co Vs are able to cleave 
the complex of ORFla and ORFlb (replicase) into 16 
mature nonstructural proteins (nsps) (Neuman et al, 
2008). Our analysis of the nsps in ORFlab revealed that 
all the three bat CoV genomic sequences contain 16 nsps 
(nspl-nspl6) in ORFlab, but the cleavage sites are dif¬ 
ferent for nsp3 or nsp5 in different Co Vs. The length of 
deduced amino acids of putative nsps, their first-last 
residue and position in replicase are shown in Supple¬ 
mentary Table S2. 

Base on the nearly complete genomic sequence ob¬ 
tained, JTAC2 possesses the same genome structure as 
PEDVs with 7 genes in the order: 5'-ORFla, lb, spike 
(S), 3a, envelope ( E ), membrane (M) and nudeocapsid 
(AO-3' (Figure 2A). JTAC2 showed the nearest relation¬ 


ship (87.9% in ORFla and 92.8% in ORFlb ) with Lushi 
MI bat CoV isolates Neixiang-14 and Neixiang-52, but 
the latter two have very limited sequences available for 
further analysis. The recent PEDV-1C isolated from a 
piglet with diarrhea and vomiting (Sun et al, 2015) was 
therefore used for sequence comparison and genomic or¬ 
ganization analysis, since it has been fully sequenced and 
shares high identity with JTAC2 ( figure 1A and Supple¬ 
mentary Table S3). The aa identity comparison shown in 
Supplementary Table S3 suggests that JTAC2 is a novel 
aCoV. 

JTMC15 is a SARSr-BatCoV having the same gen¬ 
ome organization as other SARSr-BatCoVs (e.g., Rfl), 
but sequence deletions were observed in ORFla and N, 
and between genes lb-8. A 579-nt deletion in ORFla of 
JTMC15 was also observed in SARSr-BatCoV Rs672 
from a Rhinolophus sinicus bat (Yuan et al, 2010) and a 
human SARS-CoV ShanghaiQXC2 from the late phase 
of the 2003 epidemic (GenBank #AY463060). This 579- 
nt deletion results in a 193-aa deletion of nsp3 in ORFla, 
from residues 1059 to 1251 in the nucleic acid-binding 
(NAB) domain (Serrano et al, 2009). A second deletion 
in the N gene of JTMC15 (1156-1158 nt, one residue 
Q j68 ) was also found in 3 SARSr-BatCoV strains, 
Rp/Shaanxi2011 (Yang et al, 2013), Rml (Li et al, 2005) 
and 279/2005 (Tang et al, 2006). Interestingly, four dis¬ 
continuous deletions were identified in JTMC15 between 
genes 7b and 8 , which is unique in JTMC15, resulting in 
an ORF shift and elimination of gene 8 ( figure 2B). 
Similar to known Co Vs, extensive S gene variations 
were also observed in JTMC15, resulting in low aa iden¬ 
tities with other SARSr-BatCoV strains (the highest be¬ 
ing 86.1% to Rfl) as compared with other gene frag¬ 
ments in the genome. Receptor-binding motif (RBM) is 
an extended loop that lies on the surface of the receptor 
binding domain (RBD) of the spike protein, and is the 
most important domain for SARSr-BatCoV to recognize 
its host receptor, angiotensin-converting enzyme 2 
(ACE2) (Ren et al, 2008; Baez-Santos et al, 2015). Fur¬ 
ther alignment of the deduced amino acid sequences of 
RBM (55 aa) showed a closer relationship of JTMC15 to 
SARSr-BatCoVs than to human or civet SARS-CoVs 
(Supplementary Figure SI). Taking the above altogether, 
as shown in figure 2A and Fable 2, there are 13 genes 
predicted in JTMC15: 5'-ORFla, lb, S, 3a, 3b, E, M, 6, 
7a, 7b, N, 9a, 9b-3'. Apart from gene 7b (83.0%) and S 
(86.1%) all ORFs of JTMC15 had high aa identities to 
Rfl, ranging from 94.4% (9b gene) to 99.1% (Mgene), 
indicating that JTMC15 is a new variant within the SAR¬ 
Sr-BatCoV Rfl species. 

For JPDB144, the genome organization is almost the 
same as HKU4-4, with 10 genes in the order: 5'-ORFla, 
lb, S, 3a, 3b, 3c, 3d, E, M, N- 3' ( figure 2A). However, 
two differences were observed in the nsp2 of JPDB144: 
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Figure 2. Predicted genome organizations of JTAC2, JTMC15 and JPDB144. (A) Nonstructural proteins are represen¬ 
ted by open boxes, structural proteins by filled boxes. Apostrophes in JTAC2 identify unsequenced regions. (B) Se¬ 
quence comparison showing the ORF shift of gene 7b of JTMC15 caused by the discontinuous deletions (represented 
by dots), resulting in elimination of gene 8 as compared to other SARS- and SARSr-CoVs. Nucleotide position are de¬ 
termined referencing strain BJ01. Stop codons and start codons are in bold fonts. Hu: human SARS-CoV; Ci: civet 
SARS-CoV; Bt: Bat SARS-CoV. 


a 12-nt insertion (residues 1143 to 1146 of la) and a 3-nt 
deletion (residue 1155 of la). Other JPDB144 ORFs 
were the same as HKU4-4 in length, sharing aa identit¬ 
ies of between 88.8% (3c gene) and 98.8% (E gene); 
however, an aa sequence comparison of JPDB144 ORFs 
to those of HKU5 and MERS-CoV strains in the 
Betacoronavirus lineage 3 showed rather low similarit¬ 
ies (Supplementary Table S4). 

DISCUSSION 

As shown in 7 igure 1, diverse aCoVs and PCoVs have 
been identified in the present study from different bats 
sampled at 25 locations in 4 provinces and the Tibet 
Autonomous Region, demonstrating the wide distribu¬ 
tion of CoVs among a range of bat species. Of 8 aCoVs 
identified, YDB5C is the first bat-borne CoV identified 
in the Tibet and Himalayan area, detected in 1 of 15 Hip- 
posideros cineraceus bats collected in Yadong county of 
Tibet, located at the southern edge of the Himalayas bor¬ 
dering on Bhutan and India. Another newly identified 
CoV, MLHJC4, was detected in Rhinolophus sinicus in 
Yunnan province, which phylogenetically clustered 


closely with YDB5C, both showing 94% nt identity to 
HKU2/GD/430/2006 identified in Guangdong (Lau et al, 
2007), indicating that this type of aCoV has a wide range 
of bat reservoirs and geo-distribution in south-west 
China and perhaps neighboring regions. In addition, six 
other aCoV sequences JTAC2, MLHJC2, MLHJC6, 
MLHJC8, MLHJC22 and MLHJC34, found in this study 
clustered as two novel CoV groups. Although not novel, 
the PCoVs identified here showed abundance in genetic 
and geographical diversities. It is interesting to note that 
SARSr- and MERS-like CoVs were identified, particu¬ 
larly JTMC15 isolated in Jilin province - the first SAR- 
Sr-Bat CoV to be discovered in Northeast China. 

In consideration of bat species, sampling locations and 
CoVs diversities, four pathogen/host/environment situ¬ 
ations can be proposed. First, a single bat species at one 
location (even a single cave) harboring different CoV 
species (e.g., Rhinolophus sinicus collected at the same 
site in Menglian county, Yunnan province, harboring three 
CoV species: HKU2-like, SARSr- and new aCoVs). Sec¬ 
ond, a single bat species roosting at different locations 
harboring the same CoVs (e.g., Rousettus leschenaulti 
collected in Mengla county (south Yunnan) and Wand- 
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Table 2. Comparison of ORF amino acid identities of JTMC15 and other SARS- and SARSr- CoVs # 


ORF 

JTMC15 


Rfl 


Rs672 


BJ01 

Length 

Length 

% identity 

Length 

% identity 

Length 

% identity 

la 

4185 

4378 

98.0 

4190 

93.8 

4383 

93.5 

1b 

2704 

2704 

98.1 

2704 

98.1 

2704 

98.0 

S 

1236 

1241 

86.1 

1255 

81.9 

1241 

76.7 

3a 

274 

274 

98.2 

274 

92.0 

274 

86.2 

3b 

114 

114 

97.4 

114 

91.3 

114 

90.4 

E 

76 

76 

94.8 

76 

96.1 

76 

96.1 

M 

221 

221 

99.1 

221 

98.2 

221 

97.7 

6 

63 

63 

96.9 

63 

92.2 

63 

89.1 

7a 

122 

122 

98.4 

122 

93.5 

122 

91.9 

7b 

52 

44 

83.0 

44 

84.9 

44 

79.2 

8(8af 

- 

122 

- 

121 

- 

39 

- 

-(8b) 

- 

- 

- 

- 

- 

84 

- 

N 

420 

421 

98.1 

422 

96.7 

422 

96.2 

9a 

97 

97 

94.9 

98 

79.6 

98 

79.6 

9b 

70 

70 

94.4 

70 

83.1 

70 

84.5 


Note: Abbreviation and accession numbers: Rf 1, DQ412042; Rs672, FJ588686; BJ01, AY278488. £ Gene 8 in SARSr-CoVs 
is described as 8a and 8b in SARS-CoVs. 


ing county (west Yunnan) harboring the same pCoVs). 
Third, multiple bat species sampled at the same site har¬ 
boring the same CoVs (e.g., Rousettus leschenaulti and 
Megaerops kusnotei sampled at the same location in 
Wanding county harboring the same PCoVs). Different 
bat species collected at different locations may even har¬ 
bor the same CoVs: e.g., a Hipposideros cineraceus in 
Yadong, Tibet and a Rhinolophus sinicus in Menglian, 
Yunnan harbored HKU2-like viruses, and another Rhino¬ 
lophus sinicus in Menglian and a Rhinolophus fer- 
rumequinum in Tonghua county, Jilin province, har¬ 
bored SARSr CoVs. Altogether, the data provided fur¬ 
ther evidence for the wide distribution of CoVs among 
bat populations in China, and for the suggestion that dif¬ 
ferent CoVs employ different bat species as reservoirs. 

The present study has identified genetically diverse 
bat-borne CoVs, which were detected from intestinal tis¬ 
sue specimens of different bat speices of wide geograph¬ 
ic distribution. But we failed to detect any CoV se¬ 
quence from the respiratory specimens that probably due 
to the low virus load in the lung or specific intestine- 
tropism of CoV in bats. Bats are considered the gene 
source of Alphacoronavirus and Betacoronavirus (Woo 
et al, 2012), especially of pathogenic CoVs that cause 
public threats. In last decade, increasing number of CoVs 
have been identified in bats, in which the viral genes 
were presumably originated and evolved with high muta¬ 
tion and recombination rates (Woo et al, 2007). It is ap¬ 


parent that large numbers of circulating CoVs remain 
unidentified, and are evolving worldwide within the bat 
population. Investigations in unexplored regions are 
therefore urgently needed to gain further insights into 
CoV diversity and evolutionary dynamics. 
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